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SUMMARY 


The accurate prediction of loads on flexible, low aspect-ratio aircraft is critical to the 
design of reliable and efficient vehicles. The conditions for structural design frequently 
involve nonlinear aerodynamics. 

Under previous NASA contracts (NAS1-12875, NAS1-14141, and NAS1-14962), a large 
experimental data base for three wing shapes was obtained and linear theoretical 
methods were evaluated. The current contract, NAS1-15678, extends the evaluation of 
state-of-the-art theoretical predictive methods to two separated-flow computer programs 
and also evaluates a semi-empirical method for incorporating the experimentally 
measured separated-flow effects into a linear aeroelastic analysis. 

The resultant three tasks have been documented separately. This volume describes the 
development and evaluation of a semi-empirical method to predict pressure 
distributions on a deformed wing by using an experimental data base. The evaluation of 
R. P. White’s (RASA Division of Systems Research Laboratories) separated-flow method 
(Task I) is described in NASA CR-3640. The evaluation of The Boeing Company’s 
Three-Dimensional Leading-Edge Vortex (LEV) code (Task III) is presented in NASA 
CR-3642. 

This semi-empirical method to predict pressure distributions uses an experimental data 
base in addition to a linear potential -flow solution. The experimental data accounts for 
the effects of aeroelasticity by relating the pressures to a parameter that is influenced 
by the deflected shape. Several parameters were examined before the net leading-edge 
suction coefficient was selected as the best. 



INTRODUCTION 


Accurate analytical techniques for the prediction of the magnitude and distribution of 
aeroelastic loads are required in order to achieve an optimum design of the structure of 
large flexible aircraft. Uncertainties in the characteristics of loads may result in an 
improper accounting for aeroelastic effects, leading to understrength or overweight 
designs and unacceptable fatigue life. In addition, the correct prediction of load 
distribution and the resultant structural deformation is essential to the determination 
of the aircraft stability and control characteristics, control power requirements, and 
flutter boundaries. The alternative to using satisfactory analytical techniques is the 
increased use of expensive, time-consuming wind tunnel testing for each aircraft 
configuration. 

The problem of accurate load prediction becomes particularly acute for aircraft with low 
aspect-ratio wings where critical design conditions occur in the transonic speed regime. 
In this region at typical design angles of attack, the flow is generally nonlinear - mixed 
flow, embedded shocks, separation, and vortex flow. 

A program was started in 1974 to systematically obtain experimental pressure data for 
an arrow wing throughout the subsonic, transonic, and low supersonic Mach numbers. 
This program was comprised of three NASA contracts: NAS1-12875, NASl-14141, and 
NAS1-14962 (documented in refs. 1 through 12). As the specific objective was to 
understand the change in load with aeroelastic deformation, three wing shapes were 
tested - all with the same planform and thickness distribution. The first wing was flat 
(no camber or twist); the second has a spanwise twist (typical of aeroelastic 
deformation) but no camber; and the third has the same twist with camber 
superimposed. 

In addition to the creation of a data base, which is useful for evaluating aeroelastic 
effects, a second objective was to evaluate state-of-the-art theoretical methods that 
might be used for this purpose. Primarily these methods were linear and the 
evaluations showed that linear theories are adequate at low angles of attack, which are 
typical of cruise conditions, and are basically capable of predicting loading changes due 
to smooth changes in wing shape at these low angles. However, at the higher angles of 
attack, which are typical of structural design conditions, these methods are not useful 
because the flow is nonlinear due to leading-edge separation of the flow. The limited 
comparisons that were made with advanced separated-flow methods indicated some 
hope, even though the aerodynamic panel model available at that time was very crude 
(only a few panels to represent the camber surface). 

The current evaluation of methods for predicting pressure distributions when the flow is 
separated is divided into three tasks. Two currently available computer codes were 
evaluated in Tasks I and III, and an approach involving semi-empirical corrections to 
linear theory was investigated in Task II. The three tasks are essentially independent 
efforts and are documented separately: Task II, the development and evaluation of a 
semi-empirical method in this document; Task I, an evaluation of R. P. White’s 
computer code in NASA CR-3640; and Task III, an evaluation of Boeing’s 
Three-Dimensional Leading-Edge Vortex computer code in NASA CR-3642. 

The semi-empirical methods evaluated in this task use wind tunnel pressure data from 
a representative model to predict the nonlinear aeroelastic load distribution on a wing 
when the flow is dominated by leading-edge vortices and areas of separated flow. 



Several possible parameters were examined with which to relate the net pressures. The 
most successful of these is the net leading-edge suction coefficient, which is calculated 
from the theoretical pressures and directly accounts for the effects of aeroelastic wing 
camber and twist. This is because the pressures from the appropriate shape are used to 
calculate the suction coefficient. The results were best when using the leading-edge 
suction at the same body station. 
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SYMBOLS 


b 

BL 

c 

c, M.A.C. 
C B 

C c 

C c 


'm.25c 


C N 

cl 


's,net 


D 

M 

MS 


Pt 

q 

s 


wing span, cm 

buttock line, cm; distance outboard from model plane of symmetry 

section chord length, cm 

mean aerodynamic chord length, cm 

surface bending moment coefficient referenced to y ref , positive wingtip up 
surface chord force coefficient; positive aft 
section chord force coefficient; positive aft 

surface pitching moment coefficient, referenced to 0.25 M.A.C. ; positive 
leading edge up 

section pitching moment coefficient referenced to section leading edge or 
side of body; positive leading edge up 

section pitching moment coefficient referenced to section 0.25c; positive 
leading edge up 

surface normal force coefficient; positive up 

section normal force coefficient; positive up 

. , measured pressure - reference pressure 

pressure coefficient = — 

net leading-edge suction coefficient 
body diameter, cm 
Mach number 

model station, cm; measured aft along the body centerline from the nose 
static pressure, kN/m 2 
total pressure, kN/m 2 
dynamic pressure, kN/m 2 

reference area used for surface coefficients, cm 2 
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S h area of streamwise strip associated with a pressure station, cm 2 ; used in 

summation of section force coefficients (app. B) 

x,y,z general coordinates for distances in the longitudinal, lateral, and vertical 

directions respectively 

y re f distance outboard of model centerline of the bending moment reference 

point, cm 

a corrected angle of attack, degrees; the angle between the wing-root chord 

and the relative wind measured in the model plane of symmetry; includes 
compensation for sting deflection, tunnel-flow angularities, and wall 
effects; positive nose up with respect to relative wind 

a sec wing twist angle relative to wing reference plane, degrees; positive 

leading edge up 

/3 angle of sideslip, degrees; positive nose left with respect to relative wind 

AC p NL nonlinear incremental pressure coefficient 

8 control surface deflection, degrees; positive leading edge down for leading 

edge (see exception in app. B) and trailing edge down for trailing edge 

r) fraction of wing semispan, y/(b/2) 

A sweep angle, degrees; measured from a line perpendicular to the model 

centerline, positive aft 

4> angle defining location of pressure orifices on the surface of the 

cylindrical body at a constant MS, degrees; measured from the top of the 
body 


Subscripts: 

exp. experimental data 

L.E. leading-edge control surface 

lin. linear theory 

r wing root 

s referenced to segment of local chord 

T.E. trailing-edge control surface 
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DATA BASE 


The data obtained, both experimental and theoretical, have been presented in several 
papers (refs. 1 through 3) and are presented in more detail in numerous NASA reports 
(refs. 4 through 12). 


WIND TUNNEL MODELS 

The configuration chosen for this study was a thin, low aspect-ratio, highly swept wing 
mounted below the centerline of a high fineness-ratio body. The general arrangement 
and characteristics of the model are shown in figure 1. Two complete wings were 
constructed for contract NAS1-12875, one with no camber or twist and one with no 
camber but with a spanwise twist variation. A third wing with camber and twist was 
constructed for contract NAS1-14962. Deflectable control surfaces were available on all 
three of these wings. 

The three wings, body, and fin used to create this data base are described in detail in 
appendix A. The wings all have the same planform, thickness distribution, and 
placement of orifices. The twisted wing and the cambered-twisted wing have the same 
twist, i.e., the coordinates of the leading edges and trailing edges of the two wings are 
the same. This twist distribution is shown in figure 2. Sections at the root, midspan, 
and tip of the cambered-twisted wing (fig. 3) show not only the camber but the position 
of the sections of the cambered-twisted wing and the twisted wing, relative to the wing 
reference plane (flat wing). The flat wing had a sharp-leading-edge segment in addition 
to the rounded-leading-edge segment common to all three wings. 

The capability to measure the detailed load distribution on the wing and body of this 
configuration was provided by distributing 300 pressure orifices on the model. Each 
wing had 217 pressure orifices equally divided into seven streamwise sections on the 
left half. Orifices were located on both the top and bottom surfaces at the chordwise 
locations shown in figure 4. Pressure orifices were located on the body in five 
streamwise rows of 15 orifices each. An additional eight orifices in the area of the 
wing-body junction made a total of 83 orifices on the left side of the body. 

WIND TUNNEL TESTING 

The experimental data used in this study were obtained in the Boeing Transonic Wind 
Tunnel (BTWT) under NASA contracts NAS1-12875 and NASl-14962. A description of 
the tunnel and tests are in appendix A. The current study was limited to the wings that 
had both leading-edge and trailing-edge control surfaces undeflected. Table 1 shows a 
summary of these data. 

DATA 

The measured pressures were edited, as necessary, to account for plugged or leaking 
orifices or missing data points. The pressure coefficients were then integrated, as 
described in appendix B, to obtain streamwise section coefficients and total surface 
coefficients. When pressure coefficients were required at points other than where 
measured, a linear interpolation was used. 
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Table 1. -Summary of Subsonic /Transonic Test Conditions by Test and Run Number 







































BACKGROUND 


The objective of Task II is to evaluate semi-empirical techniques for predicting pressure 
distributions in the presence of separated flows for application to the nonlinear 
aeroelastic loads analysis of highly swept wings. 

The study is specifically directed at identifying the best approach for using wind tunnel 
pressure data from a representative model to predict the load distribution on an elastic 
airplane. This loads prediction technique will be used to account for the effects of 
aeroelasticity on the airload distribution, particularly when the flow over the wing is 
dominated by leading-edge vortices and areas of separated flow. 

Wind tunnel pressure measurements have shown that the airload distributions on the 
surface of highly swept, low aspect-ratio wings at moderate to high angles of attack are 
highly nonlinear due to flow separation and the formation of leading-edge vortices. 
While a test provides accurate measurement of the airloads on a rigid model, including 
nonlinearities, there are no satisfactory methods for predicting these airloads including 
the effect of aeroelastic twist and camber. Assuming potential flow, only the linear 
loads and their linear variations due to aeroelasticity can be predicted by linear 
theories. Yet, the structural engineer who is faced with showing compliance of the 
structure to the maneuver and gust criteria, needs accurate aeroelastic load predictions 
throughout the range of airplane angle of attack. 

In 1979, Boeing Independent Research and Development money supplemented NASA 
funds for investigating semi-empirical techniques for the predicting of nonlinear wing 
pressures on a deformed wing, using a known experimental distribution on a rigid wing. 
This research has been made possible by the availability of a large, wind tunnel data 
base, which was generated under previous NASA contracts. Data are included for three 
arrow wings that were deformed to simulate aeroelastic deformation, one flat, one 
twisted, and one cambered and twisted (see figs. 1 through 4), but having identical 
planform and thickness distribution. The nonlinearities that occur in total wing lift and 
local pressures on this configuration are documented in references 4 through 12. 
Comparisons between experimental data and data predicted using several linear 
theoretical methods are shown in references 4, 7, 8, 9, and 12. 

Typical wing upper surface isobars are shown in figure 5 for the flat wing with rounded- 
leading edge at Mach = 0.85 for two angles of attack. At a = 2.1°, the flow over the 
upper surface is attached, and the pressure distribution can be predicted using linear 
theory. These theories will also accurately predict the changes in loading due to 
moderate aeroelastic deformation at these angles of attack (ref. 12). At a = 11.9°, the 
flow has separated along the leading edge, and a strong vortex dominates the flow over 
the wing. Total wing aerodynamic coefficients, wing section data, and local pressures all 
show highly nonlinear characteristics as a function of angle of attack. The linear 
theories cannot predict how such a lift distribution will vary when the wing deforms 
under load. 

Figure 6 shows the section normal force coefficient versus angle of attack for a 
streamwise wing section located at 80-percent semispan for both the flat wing and the 
cambered-twisted wing. The objective of Task II will have been met when a 
semi-empirical technique, using the data from the flat wing as a data base, accurately 
predicts the cambered-twisted wing test data. 
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Two typical methods used in the past for aeroelastic load prediction are outlined in 
figure 7. Part (a) shows the prediction of section C n for the cambered-twisted wing using 
the linear potential-flow theory FLEXSTAB; the shaded area shows the error between 
this linear prediction and actual test data. 

Figure 7 (b) shows how the flat wing test data, corrected by linear theory for the 
incremental effect of twist and camber, would compare with actual test data. Clearly, 
neither of these methods is satisfactory. 


9 



APPROACH 


The goal in this task was to find some parameter that would relate the net pressure 
coefficients from a data base wing shape to the net pressure coefficients on a deformed 
wing shape. 


PARAMETERS SELECTED FOR EVALUATION 

The following parameters were selected for examination: 

o Local section angle of attack (wing angle of attack plus local twist angle) 
o Local section normal force coefficient 
o Full theoretical leading-edge suction coefficient 

o Net leading-edge suction coefficient (includes effect of leading-edge shape and wing 
deformation) 


For the first three of these parameters, the full net pressure coefficient would be 
obtained directly. However, only the nonlinear portion of the net pressure coefficient 
would be obtained from the data base when using the net leading-edge suction 
coefficient, and it would be added to the net pressure coefficient calculated using linear 
theory. 

The use of the leading-edge suction analogy of Polhamus for predicting the gross effects 
of leading-edge vortices has been well substantiated for very sharp, thin, plate-type 
wind tunnel model wings (see ref. 13). This analogy uses the full theoretical 
leading-edge suction coefficient, which, as described in appendix C, is obtained from the 
net pressure distribution predicted by FLEXSTAB (refs. 14 through 17). 

Experimental data have shown that wing thickness affects the growth and development 
of leading-edge vortices. Furthermore, the effects of wing thickness on the leading-edge 
vortex development depend on whether the leading edge is sharp or rounded. 

The concept of net leading-edge suction has been formulated to account for thickness 
and leading-edge shape effects on leading-edge vortices (see refs. 18 and 19). The net 
suction coefficient relates the local vortex strength to the leading-edge suction 
calculated at what amounts to a reduced (or effective) angle of attack. This reduced 
angle of attack for wings with sharp-nose airfoils is equal to the average angle of attack 
over the front portion of the local streamwise airfoil upper surface. The vortex growth is 
assumed to begin when this angle becomes positive. 

Using the same terminology, the effective angle of attack of airfoils with rounded noses 
is also reduced. At each section, the effective angle of attack is defined as the difference 
between the local angle of attack and the angle at which the leading-edge suction 
equals the parabolic nose drag (the local nose is approximated by a parabola). It is 
assumed that the local vortex growth begins when this angle becomes positive. The 
procedure used for this calculation is developed in references 18 and 19 and is 
summarized in appendix C. 
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The relationship of net leading-edge suction coefficient C 8 net to angle of attack is shown 
in figure 8 for the flat wing, the twisted wing, and the cambered- twisted wing. The 
square root of C 8 net is used because it is a linear function of the angle of attack in the 
range shown (-8° to +16°). 


EVALUATION CRITERIA 

For each method examined under this task, the following points were considered: 
o Definition of a theoretical basis for the method 
o Identification of the key assumptions 

o Determination of how the method would be used with test data in an aeroelastic 
analysis 

o Evaluation of the method by comparing these predictions to linear-theory 
predictions and to test data for the appropriate configuration from references 1 
through 12 

- Predict pressure distributions on the cambered-twisted wing using the flat wing 
data for the data base 

- Predict pressure distributions on the cambered-twisted wing using the twisted 
wing data for the data base 

- Predict pressure distributions on the twisted wing using the flat wing data for 
the data base 

o Identification of method limitations, deficiencies, improvement options, and/or 
alternate methods. 
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EVALUATION OF PARAMETERS 


All evaluations were made at a Mach number of 0.85, as the data at that Mach number 
appears to be representative of the data on this configuration. 

LOCAL SECTION ANGLE OF ATTACK 

Relating pressure distributions to the local section angle of attack clearly could not 
account for the effects of wing camber changes. Therefore, the examination of this 
approach was limited to plotting of some of the data. 

SECTION NORMAL FORCE 


KEY ASSUMPTION 

In this method, it must be assumed that aeroelastic wing twist and camber changes will 
not appreciably affect the shape of the chordwise pressure distributions for a specific 
section normal force coefficient and spanwise station. 

APPLICABILITY FOR USE IN AN AEROELASTIC ANALYSIS 

An aeroelastic method, based on relating the local net pressure to the section normal 
force coefficient, could be implemented as shown in figure 9. This method would be used 
to predict aeroelastic twist effects by using wind tunnel model data in the following 
manner: 

1. Establish a data base using experimental pressure data. At each spanwise station, 
the net pressure coefficients are related to the section normal force coefficient 
obtained by the integration of these pressures. 

2. Set external pressure field for FLEXSTAB to zero. 

3. Execute FLEXSTAB - a trimmed, flexible airplane solution - for the required load 
factor. Include the appropriate external pressure field. 

o Evaluate theoretical pressure distributions for this deformed shape 

o Integrate the pressure coefficients (include external pressure field) to obtain 
the section normal force coefficient 

4. Obtain the net pressure distribution from the data base for the section normal 
force coefficients, calculated in step 3, to use as the external pressure field in the 
FLEXSTAB solution. 

5. Repeat steps 3 and 4 until the following converge: 
o airplane angle of attack and, 

o elevator angle required to trim 

As stated above, this approach is based on the assumption that the shape of the 
chordwise pressure distribution is not affected by aeroelastic deflections. Since 
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chordwise camber does affect the shape of the pressure distribution, this approach is not 
expected to be applicable when there are chordwise camber changes. 

EVALUATION 

In order to evaluate this parameter, comparisons were made at a number of spanwise 
stations of the net pressure distributions that produced the same section normal force 
coefficients for the flat and twisted wings. This comparison is shown in figure 10. A 
similar comparison of data on the twisted and cambered-twisted wings is shown in 
figure 11. 

The shapes of these pressure distributions on the flat wing and the twisted wing show 
small differences that would affect the ability to predict local loading and load paths. 
The pressures on the twisted wing and the cambered-twisted wing indicate that the 
simulated camber changes affect the shape of the chordwise load distributions 
significantly over most of the wing and, therefore, would affect not only local loading 
but overall torsion. In order to obtain a partial assessment of the effects of the 
difference in chordwise pressure distribution on the overall load distribution for the flat 
and the twisted wings, calculations were made of the bending moment, shear, and 
torsion about the reference axis; this was done by integrating the pressures shown in 
figure 10. 

The results indicate that the difference in the pressure distribution, due to twist, have a 
very small effect on the overall wing shear, bending moment, and torsion. However, this 
can be misleading when analyzing the structure of low aspect-ratio wings that generally 
have multiple load paths. Similar calculations for the camber effect clearly show a 
significant discrepancy in torsion, which in turn would affect the aeroelastic behavior of 
the wing. 

METHOD LIMITATIONS 

In addition to the inaccuracies mentioned above, two major difficulties have been 
identified that affect the validity and potential usefulness of using normal force as a 
means of predicting pressure distributions. 

1. This method requires the calculation of changes in section normal force caused by 
aeroelastic distortion. The change in aeroelastic shape affects both the strength 
and location of the leading-edge vortices. Hence, the effect of those changes on the 
spanwise normal force distribution is highly nonlinear and cannot be predicted by 
a linear-theory analysis, such as FLEXSTAB. 

2. A second and perhaps more serious difficulty is illustrated by the section normal 
force comparisons in figure 12. At the inboard stations, the section normal force 
coefficients for the study wings increased with angle of attack over the range 
plotted. However, near the wing tip, the section normal force coefficients initially 
increase as the leading-edge vortex develops; then as the leading-edge vortex grows 
and the vortex core moves inboard, the normal force drops and remains nearly 
constant. The same normal force coefficient may occur at two or more angles of 
attack; however, the chordwise load distributions usually are significantly 
different. Also, comparisons of section normal force indicate the possibility of 
having a particular section normal force within the linear-flow region on the flat 
wing, and within a highly nonlinear-flow region on the cambered-twisted wing. 
Again, the chordwise load distributions are significantly different. 
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FULL THEORETICAL LEADING-EDGE SUCTION COEFFICIENT 

The possibility of relating net pressures to the full theoretical leading-edge suction was 
briefly explored. The concept of using the net leading-edge suction as a parameter for 
predicting the appearance, growth, strength, and location of the leading-edge vortex, 
and hence, its incremental, nonlinear effect on the pressure distribution, seemed to be 
more in keeping with the Polhamus suction analogy. 

NET LEADING-EDGE SUCTION COEFFICIENT 

A two-part examination was undertaken to see if a possible relationship existed 
between the net leading-edge suction coefficient and the incremental net pressures. The 
first part of the examination used the more obvious approach, which was to relate these 
pressures to the leading-edge suction at the same spanwise location, with sections as 
defined in figure 13. The second approach (considering the flow on such a highly swept 
wing is frequently described as conical) was to relate the pressures to leading-edge 
suction at the same body station as shown in figure 14. In both cases, the leading-edge 
suction coefficient was evaluated using pressures in the streamwise sections that are 
typical of wing paneling for FLEXSTAB. 

KEY ASSUMPTION 

This empirical method assumes that the net incremental nonlinear pressures are a 
function of the net leading-edge suction coefficient C s net as calculated from the 
theoretical FLEXSTAB pressures. This would imply that the normal force and pitching 
moment coefficients resulting from these pressure distributions are also a unique 
function of C s net . These assumptions hold for either the streamwise or constant body 
station method. 

APPLICABILITY FOR USE IN AN AEROELASTIC ANALYSIS 

An aeroelastic method based on this type of relationship would be implemented as 
shown in figure 15. Aeroelastic loads would be predicted in the following manner: 

1. Establish a data base, as shown in figure 16, using experimental pressure data and 
FLEXSTAB-generated theoretical pressures for the same model. 

o Obtain incremental, nonlinear net pressure coefficients (difference between 
experiment and linear theory) as a function of angle of attack. 

o Calculate (see app. C) net leading-edge suction coefficients as a function of 
spanwise location and angle of attack. 

o Relate incremental, nonlinear net pressure coefficients to net leading-edge 
suction coefficient at the same angle of attack. 

2. Set external pressure field for FLEXSTAB to zero. 
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3. Execute FLEXSTAB - a trimmed, flexible airplane solution - for the required load 
factor. Include the appropriate external pressure field. 

o Evaluate theoretical pressure distributions for this deformed shape. 

o Calculate net leading-edge suction coefficients using these pressures. 

4. Obtain the incremental, nonlinear net pressure distribution from the data base for 
the net leading-edge suction coefficients calculated in step 3 to use as the external 
pressure field in the FLEXSTAB solution. 

5. Repeat steps 3 and 4 until the following converge: 
o airplane angle of attack and, 

o elevator angle required to trim. 


EVALUATION 
Streamwise Sections 

The nonlinear portions of the local net pressures for the flat, the twisted, and the 
cambered-twisted wings were displayed versus the square root of the net leading-edge 
suction coefficient at several streamwise sections. The plots tended to coincide, 
especially for the inboard sections in the vicinity of the leading edge; but no measure of 
goodness in this relationship could be shown using that particular mode. 

Section C n and C m 25c for the cambered-twisted wing were predicted, using the data 
base generated from the flat wing data, in the following manner: 

o Integrate the FLEXSTAB-generated linear pressure distribution to obtain section 
C n and C m 25c for the cambered-twisted wing. 

o Calculate the net leading-edge suction coefficient using these pressure 
distributions. 

o Obtain the incremental, nonlinear section C n and C m2 5 C from the data base for 
these suction coefficients. 

o Combine the coefficients from linear theory and from the data base. 

The predicted section C n and C m 25c were then compared to the actual test values at the 
seven streamwise sections along the span at which pressures were measured on the test 
models. This comparison of the section coefficients as a function of C 8 net is shown in 
figure 17. These coefficients are also shown as a function of angle of attack in figure 18, 
with those predicted by FLEXSTAB added to show the results of a typical linear theory. 

As shown in figure 19, the procedure was repeated for the twisted wing - again using 
the data base for the flat wing in order to isolate the twist effect. To isolate the camber 
effect, the cambered-twisted wing section coefficients were predicted by using the data 
base for the twisted wing. The results are shown in figure 20. 
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Significant errors between predicted and actual section C n and C m 2 5 C occur on the 
outboard half of the wing. Such results are clearly unsatisfactory. 

Crossflow Sections 

The evaluation procedure for crossflow sections was similar to that for streamwise 
sections. In this case, however, the net pressures were related to the leading-edge 
suction coefficients at the same body station (crossflow section). Plots of C p nl versus 
C 8 net were made for a matrix of stations and x/c locations for the three wings using the 
procedure shown in figure 16. Typical results for the three wings are shown in 
figure 21. 

Experimental normal force and moment coefficients versus angle of attack for the flat, 
the twisted, and the cambered-twisted arrow wings are shown in figure 22 for a typical 
wing section that is perpendicular to the model plane of symmetry. The normal force 
coefficient is separated into linear and nonlinear parts; it was then plotted as a function 
of the square root of the predicted net leading-edge suction coefficient for that wing. In 
this angle-of-attack range, the leading-edge suction coefficient is essentially linear with 
the square of a . 

We then considered the problem of predicting wing loads for the aeroelastically 
distorted wing. This would have to be done using experimental data for a base wing of 
the same planform, thickness distribution, and leading-edge radius. The linear portion 
of the load could be predicted for the distorted wing using a theoretical method such as 
FLEXSTAB. The nonlinear portion of the load would be obtained for the distorted wing 
from the base wing data at the appropriate leading-edge suction coefficients. The sum of 
these two parts is the predicted load. Figures 23 through 25 show, for nine crossflow 
sections, a comparison of the predicted values versus the actual test values of section C n 
and C m . Respectively, the figures show the ability to predict the effects of wing twist, 
camber, and combined twist and camber. The results shown in these figures are better 
than those shown in figures 19, 20 and 17, where the evaluation of the method using 
streamwise sections was shown. 

Figure 26 shows the comparison of predicted section coefficients for the 
cambered-twisted wing to test data and theoretical predictions as a function of angle of 
attack. These predicted section coefficients are based on the data base for the flat wing. 
At an angle of attack of 12°, the associated pressure distribution at a streamwise, 
midspan section is shown in figure 27. Again, this approach shows very promising 
results. 

The crossflow section correlation of local lift with theoretical leading-edge suction is a 
promising method for predicting the effects of aeroelastic twist and camber for an arrow 
wing that is in the presence of separated flow and a leading-edge vortex. 
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CONCLUSIONS 


PRESSURES RELATED TO SECTION NORMAL FORCE 

The net chordwise pressures for the flat wing and the twisted wing were virtually 
identical when compared at the same section normal force. This suggested the use of 
section normal force as a parameter to account for the effect of aeroelastic twist. 
However, the comparison of the twisted wing pressures and the cambered-twisted wing 
pressures was not nearly as good. This could be expected, since an airfoil camber change 
does affect the shape of the chordwise pressure distribution and as a result, affects wing 
torsion and aeroelastic deformation. Also, the pressure distributions were not unique 
with section C n . These and other limitations led to abandoning this method in favor of 
using the leading-edge suction coefficient. 

NONLINEAR PRESSURE INCREMENTS RELATED TO NET 
LEADING-EDGE SUCTION 

Work performed in this contract (see refs. 18 and 19) expanded the Polhamus 
leading-edge suction analogy derived for sharp leading-edge wings to include the effects 
of leading-edge radius and of airfoil thickness and camber on the development and 
growth of leading- edge vortices. 

STREAMWISE SECTIONS 

It appeared that linking the lift nonlinearities to the local net leading-edge suction 
coefficient, which was calculated using FLEXSTAB theoretical pressures (see refs. 14 
through 17), was the most logical and promising approach. As suggested by Polhamus, 
the pressure nonlinearities at a given streamwise wing section were related to the net 
suction at the leading edge of that same section. The loads predicted for the 
cambered-twisted wing, using a data base created from flat wing pressure data, did not 
satisfactorily match the actual loads for the cambered-twisted wing. Since the flow 
around an arrow wing is conical with respect to the wing apex, it appeared that relating 
the loads to the suction at crossflow (constant body station) sections, rather than in 
streamwise sections, would work better. 

CROSSFLOW SECTIONS 

Using this approach, the predicted pressures and resulting crossflow section lift and 
moments accounted well for the effects of twist change, camber change, and for the 
combined effect of both twist and camber when compared to measured data. These 
results represented a quantum improvement over pure theoretical predictions. The use 
of suction coefficients at crossflow sections was, therefore, selected as the best approach 
for using wind tunnel pressure data from a representative model to predict the load 
distribution on an elastic airplane in the presence of separated flows. 
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RECOMMENDATIONS 


Further substantiation of the method using net leading-edge suction at crossflow 
sections for the arrow wing should be carried out. Test data are available for these 
evaluations at Mach numbers of 0.40 to 2.50. 

The effect of moderate leading-edge deflection, regarded as local camber change, was 
thought to be predictable using this method; but a closer look at the discontinuities 
introduced at the hinge line in the theoretical FLEXSTAB pressure distribution, makes 
this approach doubtful. However, a continuous-curvature leading-edge deflection, like 
that of the mission-adaptive wing of the AFTI-111 airplane, should be a good candidate 
for further study. 

The applicability of the method to other wing configurations is of great interest, and the 
AFTI-111 variable sweep wing could provide more insight in the effects of: 

- wing sweep angle 

cranked leading edge 

The use of far-field theory to improve the leading-edge suction calculation should be 
examined. The procedure would be to scale the local leading-edge suction as calculated 
using the FLEXSTAB pressures, so that its summation along the leading edge matches 
the global suction predicted by far-field theory. 

Boeing Commercial Airplane Company 
P.O. Box 3707 

Seattle, Washington 98124 
May 1982 
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APPENDIX A 


DESCRIPTION OF DATA BASE 


WIND TUNNEL MODELS 

The configuration chosen for this study was a thin, low aspect-ratio, highly swept wing 
mounted below the centerline of a high fineness-ratio body. The general arrangement 
and characteristics of the model are shown in figure A-l. Two complete wings were 
constructed for contract NAS1-12875, one with no camber or twist and one with no 
camber but with a spanwise twist variation. A third wing with camber and twist was 
constructed for contract NAS1-14962. Deflectable control surfaces were available on all 
three of these wings. 

FLAT WING 

The mean surface of the flat wing is the wing reference plane. The nondimensional wing 
thickness distributions (shown in table A-l) deviate slightly from a constant for all 
streamwise sections to satisfy a manufacturing requirement for a finite thickness of 
0.0254 cm (0.01 in.) at the trailing edge. The wing was designed with a full-span, 
25-percent chord, trailing-edge control surface. Sets of fixed angle brackets allowed 
streamwise deflections of ±4.1°, ±8.3°, ±17.7°, and ±30.2°, as well as 0.0°. A removable 
full-span leading-edge control surface (15 percent of streamwise chord) could be placed 
in an undeflected position and also drooped 5.1° and 12.8° with fixed angle brackets. 
Both the leading- and trailing-edge control surfaces extended from the side of the body 
(0.087 b/2) to the wingtip and were split near midspan (0.570 b/2). The inboard and 
outboard portions of the control surfaces were able to be deflected separately and were 
rotated about points in the wing reference plane. An additional leading-edge control 
surface for this wing was constructed with a sharp (20° included angle) leading edge to 
examine the effects of leading-edge shape. The surface ordinates and slopes of this 
leading-edge segment were continuous with those of the flat wing at the leading-edge 
hingeline (table A-l). The sharp leading edge was smoothly faired from 0.180 b/2 into 
the fixed portion of the rounded leading edge at 0.090 b/2. 

TWISTED WING 

The mean surface of the twisted wing was generated by rotating the streamwise section 
chord lines about the 75-percent local chord points (trailing-edge control surface 
hingeline). The spanwise variation of twist is shown in figure A-2. The hingeline was 
straight and located in the wing reference plane at its inboard end (0.087 b/2) and 
2.261 cm (0.890 in.) above the wing reference plane at the wingtip. The airfoil thickness 
distribution (table A-l) and the trailing-edge control surface location and available 
deflections were identical to those of the flat wing. 

CAMBERED-TWISTED WING 

The mean surface of the cambered-twisted wing was generated by superimposing a 
camber on the twisted-wing definition but keeping the coordinates of the leading edge 
and trailing edge of the cambered-twisted wing the same as those of the twisted wing. 
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Table A- l.—Wing Half-Thickness Distribution, Percent Chord 
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The camber is defined analytically in two parts: a typical cruise airfoil (basic) camber; 
and, an estimate of the aeroelastic deformation at a moderate positive angle of attack. 
The aeroelastic deformation was based on calculations - using a typical configuration - 
of deformation under load. This definition was modified slightly to provide zero camber 
at the model centerline so this wing would fit on the existing model body. This was 
achieved by using a factor k (fig. A-3) on the basic camber term, which provides a 
transition from no camber at the model centerline to the definition camber at 0.25 b/2. 
The defined wing is smoother than indicated in this figure as section geometry was 
directly calculated at only those sections marked in figure A-3. The full equation for the 
camber is: 



The resulting nondimensional camber is shown in table A-2. The camber at the tip is 
approximately a 6° arc of a circle with the leading and trailing edges up. Sections at the 
root, midspan, and tip (fig. A-3) show not only the camber but the position of the 
sections of the cambered-twisted wing and the twisted wing relative to the wing 
reference plane (flat wing). The airfoil thickness distribution (table A-l) and the 
trailing-edge control surface location and available deflections were identical to those of 
the flat wing. 

BODY 

The body was circular in cross section and had a straight centerline. The body geometry 
is shown in figure A-l. The sting was an integral part of the model body. 


RELATIVE WING AND BODY LOCATION 

The wing reference plane was located 3.149 cm (1.240 in.) below and parallel to the body 
centerline (zero incidence). The apex (extension of the wing leading edge to the 
centerline) of the wing was located 33.496 cm (13.187 in.) aft of the model nose. 
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Table A-2.—Wing Section Camber Definition, Percent Chord 
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WING FIN 


The wing fin is a 3-percent biconvex airfoil placed streamwise and perpendicular to the 
wing reference plane on the upper surface of the cambered-twisted wing at 0.725 
semispan. The dimensions of the fin and its relationship to the wing are shown in 
figure A-4. To obtain configurations with the outboard trailing-edge control surface 
deflected, the fin is extended down to touch the top of the control surface. (See fig. A-4.) 

PRESSURE ORIFICE LOCATIONS 

All pressure orifices were located on the left side of the model and distributed as shown 
in figure A-5 and tables A-3 and A-4. The flat wing with rounded leading edge, the 
twisted wing, and the cambered-twisted wing each had 214 orifices distributed in 
streamwise pressure stations of 31 (or 30) orifices at each of seven spanwise locations. 
One of these orifices was located at the leading edge; the remainder were distributed so 
that upper- and lower-surface orifices were located at the same chordwise locations. The 
orifice locations on the sharp leading edge were identical except for the omission of the 
leading-edge orifice at each spanwise station. The 83 orifices on the body were located 
at 15 stations along the length of the model. At each station, orifices were located at 
angles of 0°, 45°, 90°, 135°, and 180° measured from the top of the body. In the area of 
the wing-body intersection, the orifices that are nominally identified as being at 135° 
and 180° were located on the wing lower surface at the same lateral location as the 
orifices at 45° and 0°, respectively. Eight additional orifices were placed on the body, 
close to the juncture of the body with the wing upper surface. 

DESIGN AND CONSTRUCTION 

The objectives of this study dictated that the contours and physical characteristics of 
the flat wing, the twisted wing, and the cambered-twisted wing be as nearly identical as 
possible. The model was constructed of steel to minimize aeroelastic deflections and to 
provide strength for testing to a Mach number of 3.0. The aft body was flared 
approximately 4° from 194.310 cm (76.500 in.) aft of the nose to provide the required 
safety factor on predicted loads (fig. A-l). The model size was selected as the best 
compromise between minimizing potential tunnel blockage and providing adequate 
room to install orifices in the model. 

A computerized lofting program was used to provide the wing definition. This definition 
was then used to machine the model components using numerically controlled 
machines. The tolerance on the contour was +0.1524, -0.0 mm ( + 0.006, -0.0 in.). The 
leading- and trailing-edge control surfaces were cut from the wings after they had been 
machined to final contour. A cut along the 15-percent chord line of the twisted wing 
removed enough material to simulate the elastic characteristics of the flat wing 
(fig. A-6). As a result of the previous tests it was determined that it was not necessary 
to remove this material on the cambered-twisted wing as the wings were very rigid. 
Fixed angle brackets (arranged as shown in fig. A-6) were used to obtain the required 
control surface deflections with all pivot points located midway between the upper and 
lower surfaces at the hingelines. The brackets were also machined on numerically 
controlled machines. The same sets of trailing-edge brackets were used on all three 
wings, and the same sets of leading-edge brackets were used for both the rounded and 
sharp leading edges. 
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Table A-3.—Wing Pressure Orifice Locations, Percent Local Chord 

(a) Section at 0.09— , chord = 102.89 cm 
2 



Nominal 

Flat wing, 

a = 0.0° 
sec 

Twisted wing, 
a sec = -°.01° 

Cambered-twisted wing, 
“sec = -°.01° 

Rounded leading edge 

Sharp leading edge 

Rounded leading edge 

Rounded leading edge 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

0.00 

0.00 



_ 

0.00 


0.00 


2.50 

2.45 

2.59 

2.61 

2.54 

2.26 

2.26 

2.58 

2.51 

5.00 

4.95 

5.07 

5.06 

5.03 

4.76 

4.76 

5.10 

5.04 

8.50 

8.45 

8.53 

8.59 

8.58 

8.40 

8.26 

8.64 

8.56 

11.30 

— 

— 

- 

11.31 

— 

— 

- 

- 

12.25 

— 

— 

— 

— 

12.23 

12.27 

— 

— 

12.50 

12.45 

12.55 

12.58 

— 

— 

— 

12.63 

12.54 

17.50 

17.49 

17.62 



17.59 

17.66 

17.64 

17.55 

20.00 

19.94 

20.08 



20.03 

20.03 

20.14 

20.00 

30.00 

29.92 

30.09 



29.98 

29.89 

30.14 

30.00 

45.00 

45.00 

45.07 



44.96 

44.89 

45.12 

45.03 

60.00 

59.98 

60.08 



60.01 

59.97 

60.11 

60.00 

70.00 

70.03 

70.13 



70.05 

69.95 

70.09 

70.04 

72.50 

72.55 

72.60 



72.58 

72.51 

72.62 

72.54 

mem 

77.53 

77.62 



77.56 

77.51 

77.63 

77.52 

■is 

85.11 

85.14 



85.03 

85.00 

85.12 

85.04 

■gin 

90.10 

90.10 



90.04 

89.98 

90.12 

90.00 

95.00 

95.09 

95.05 



94.96 

94.98 

95.10 

95.03 
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Flat wing, 

>c= 0 -°° 

Rounded leading edge Sharp leading edge 


Twisted wing, 

a = _0 47° 
sec 


Rounded leading edge 


Upper 

surface 

Lower 

surface 

0.00 

2.59 

5.05 

8.54 

12.54 

2.69 

5.00 

8.59 

12.49 

17.63 

17.61 

20.08 

20.07 

30.04 

30.09 

45.08 

45.09 

60.02 

60.13 

70.11 

70.13 

72.63 

72.61 

77.59 

77.65 

85.07 

85.13 

90.14 

90.11 

95.14 

95.10 


Cambered-twisted wing. 

a sec 

-0.47° 

Rounded leading edge 

Upper 

Lower 

surface 

surface 
































Table A-3. —( Continued) 


(c) Section at 0.35— chord = 76.69 cm 
2 ' 



Nominal 

Flat wing, 
“sec = 0 -0° 

Rounded leading edge 

Sharp leading edge 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

0.00 

0.00 


_ 

_ 

2.50 

2.45 

2.59 

2.59 

2.58 

5.00 

4.93 

5.07 

5.11 

5.04 

8.50 

8.60 

8.54 

8.65 

8.63 

10.50 

— 

— 

— 

10.46 

11.00 

— 

11.03 

— 

— 

12.50 

12.37 

- 

12.57 

— 

17.50 

17.64 

17.63 



20.00 

20.00 

20.09 



30.00 

30.01 

30.10 



45.00 

44.99 

45.09 



60.00 

60.03 

60.08 



70.00 

70.07 

70.08 



72.50 

72.55 

72.58 






77.61 

85.14 

90.09 

95.09 




Twisted wing, 
tt sec=- 1 - 70 ° 

Cambered-twisted wing, 
I a sec = —1-70° 

Rounded leading edge 

| 

j Rounded leading edge 

Upper 

Lower 

Upper 

Lower 

surface 

surface 

surface 

surface 

0.00 

2.39 

2.33 

0.00 

2.76 

2.60 

5.12 

4.78 

5.05 

5.10 

8.49 

8.32 

8.68 

8.70 

— 

— 

— 

— 

— 

— 

— 

— 

12.50 

12.33 

12.59 

12.68 

17.54 

17.53 

17.64 

17.62 

19.94 

19.84 

20.03 

20.07 

29.88 

29.87 

30.00 

29.93 

44.96 

44.79 

45.00 

45.13 

59.97 

59.89 

60.00 

60.10 

70.03 

69.90 

70.04 

70.03 

72.56 

72.44 

72.61 

72.52 

77.54 

77.51 

77.50 

77.60 

85.08 

84.96 

85.09 

84.93 

89.89 

89.89 

89.98 

90.04 

94.95 

94.86 

94.98 

95.10 


95.00 































Table A -3. —( Continued) 



Twisted wing, 
“sec = - 2 - 85 ° 

Cambered-twisted wing, 
“sec = - 2 -85° 

Rounded leading edge 

Rounded leading edge 

Upper 

Lower 

Upper 

Lower 

surface 

surface 

surface 

surface 

0.00 


0.00 


2.44 

2.38 

2.78 

2.62 

4.92 

4.80 

5.13 

5.15 

8.46 

8.38 

8.64 

8.56 

12.50 

12.31 

12.71 

12.55 

17.54 

17.24 

17.71 

17.44 

19.92 

19.83 

20.15 

19.89 

29.91 

29.85 

30.04 

29.72 

45.00 

44.85 

44.95 

44.97 

59.95 

59.92 

59.96 

59.94 

70.03 

69.88 

69.93 

69.86 

72.56 

72.44 

72.53 

72.34 

77.61 

77.43 

77.58 

77.43 

84.85 

84.90 

84.96 

84.92 

89.93 

89.93 

89.94 

89.91 

94.88 

94.93 

94.98 

94.88 
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Table A-3. —( Continued) 


(e) 


Section at 0.65—, chord = 46.46 cm 
2 



Nominal 

Flat wing, 

a = 0.0° 
sec 

Twisted wing, 

a = -3.59° 
sec 

Cambered-twisted wing, 

a = -3.59° 
sec 

Rounded leading edge 

Sharp leading edge 

Rounded leading edge 


Rounded leading edge 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

0.00 

0.00 




0.00 


0.00 


2.50 

2.56 

2.66 

2.49 

2.38 

2.18 

2.49 

2.76 

2.79 

5.00 

5.06 

5.12 

4.94 

4.95 

4.76 

5.01 

5.35 

5.40 1 

8.50 

8.55 

8.55 

8.46 

8.40 

8.32 

8.45 

8.85 

8.74 

12.20 

— 

- 

12.12 

— 

12.21 

- 

- 

— 

12.60 

12.57 

— 

— 

— 

— 

— 

12.71 

12.77 

17.50 

17.60 

17.65 



17.24 

17.44 

17.74 

17.58 

20.00 

20.17 

20.11 



19.70 

19.88 

20.19 

19.96 

30.00 

30.05 

30.11 



30.26 

29.73 

30.13 

29.85 

45.00 

45.16 

45.23 



44.75 

44.89 

45.03 

44.75 

60.00 

60.13 

60.13 



59.81 

59.87 

60.02 

59.99 

70.00 

69.89 

70.12 



69.92 

69.90 

70.09 

69.88 

72.50 

72.59 

72.69 



72.38 

72.49 

72.83 

72.15 

77.50 

77.74 

77.76 



77.22 

77.49 

77.56 

77.43 

85.00 

85.25 

85.32 



84.79 

84.93 

84.93 

84.76 

90.00 

90.22 

90.21 



89.70 

89.92 

89.95 

89.98 

95.00 

95.13 

95.27 



95.12 

94.86 

94.97 

94.98 
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Nominal 


Flat 

a 

sec 

wing, 
= 0.0° 


Twisted wing, 
“sec = - 3 -84° 

Cambered-twisted wing, 
^c = -3-84° 

Rounded leading edge 

Sharp leading edge 

Rounded leading edge 

Rounded leading edge 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

0.00 

0.00 




0.00 


0.00 


2.50 

2.55 

2.47 

2.50 

2.46 

2.33 

2.43 

2.76 

2.62 

5.00 

5.01 

5.02 

5.01 

4.93 

4.86 

4.74 

5.27 

5.21 

8.50 

8.55 

8.59 

8.58 

8.41 

8.32 

— 

8.78 

8.54 

12.50 

12.50 

‘ 

12.58 


12.47 

12.43 

12.69 

12.58 

17.50 

17.53 

17.57 



17.36 

17.47 

17.83 

17.34 

20.00 

20.16 

20.13 



19.79 

19.82 

20.11 

19.79 

30.00 

30.00 

30.11 



29.83 

29.83 

30.15 

29.48 

45.00 

44.91 

45.15 



44.81 

44.91 

44.81 

44.75 

60.00 

59.94 

60.10 



59.80 

59.92 

59.84 

59.79 

70.00 

70.06 

70.11 



69.89 

69.87 

69.77 

69.94 

72.50 

72.61 

72.60 



72.22 

72.39 

72.50 

72.33 

77.50 

77.73 

77.72 



77.29 

77.41 

77.22 

77.40 

85.00 

85.25 

85.18 



84.80 

84.95 

84.92 

84.92 

90.00 

90.20 

90.34 



90.62 

90.03 

90.19 

90.09 

95.00 

95.41 

95.49 



95.71 

95.00 

95.05 

94.94 






































Table A-3.— (Concluded) 


<g> 


Section at 0.93—, chord = 18.25 cm 
2 


T.E. hingeline 



Nomina 

Flat wing, 
«sec = 0 -°° 

Twisted wing, 
“sec = - 4 - 14 ° 

Cambered-twisted wing, 
a =—4 14° 

sec 

Rounded leading edge 

Sharp leading edge 

Rounded leading edge 

Rounded leading edge 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

2.26 

4.79 

8.13 

12.16 

0.00 

2.51 

5.00 

8.50 

11.59 

12.25 

0.00 

1.70 

4.38 

7.89 

12.33 

1.81 

4.68 

8.24 

2.12 

4.72 

8.21 

12.19 

1.86 

4.52 

8.06 

— - 

0.00 

1.74 

4.41 

7.92 

11.59 

2.59 

4.65 

8.23 

0.00 

2.77 

5.11 

8.64 

12.64 

mm 

17.36 

16.60 



16.60 

17.49 

18.03 

16.83 

— 

19.78 

19.81 



19.58 

19.96 

19.94 

19.44 

— 

29.67 

29.00 



29.17 

29.62 

30.22 

28.66 

Bl 

44.70 

44.80 



44.12 

44.44 

44.33 

44.77 

all 

59.68 

59.47 



59.18 

59.71 

59.47 

59.38 

70.00 

69.69 

70.33 



68.99 

69.31 

69.10 

70.07 

72.50 

72.15 

71.89 



71.59 

72.01 

71.78 

72.74 

77.50 

77.38 

77.31 



76.80 

77.12 

76.49 

77.36 


84.62 

84.90 



84.54 

84.82 

84.93 

85.29 


89.51 

89.81 



89.21 

89.74 

90.72 

90.35 

j 95.00 

94.46 

94.68 



94.41 

94.56 

95.26 

94.87 
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Table A-4.-Body Pressure Orifice Locations 
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For the first 149 runs of the first test in the BTWT, pressure readings at these orifices did not always stabilize. 


Pressure tubing used in this model was 1.016 mm (0.040 in.) o.d. Monel with a 
0.1524 mm (0.006 in.) wall thickness. The major channels for wing pressure tubing were 
machined into the surface. The detailed grooves required to route tubing from the 
orifices to these channels were cut by hand. The pressure orifices were installed normal 
to and flush with the local surface. After installation of the pressure tubing, the grooves 
were filled with solder and brought back to contour by hand-filing to match templates 
prepared by numerically controlled machining. 

Quick disconnects were used at the wing-body junction to reduce the time required for 
installing a different wing. Unfortunately, by the time the cambered-twisted wing was 
installed in the test section, one quick-disconnect block had become worn out due to the 
two previous tests and model checkout. The connection did not seal properly and 
measurements at a series of orifices (x/c from 0.125 through 0.600) on the lower surface 
at 0.80 b/2 were not sufficiently accurate to be used. Data values to be used in the 
integration were obtained by linear spanwise interpolation between adjacent sections. 

The tubing for body pressure orifices was run through the hollow center of the model 
body rather than running it in grooves in the outside contour. Tubing from all the 
orifices was routed through the hollow body to the scanivalves located in the body nose. 
Wiring from the scanivalves was routed through the body to the sting. 


The nose portion of the body was removable to provide access to the fifteen 24-position 
scanivalves. Figure A-l shows the aft body location of the strain gages that were used to 
measure normal force and pitching moment. 

PRESSURE INSTRUMENTATION 

The model was instrumented with fifteen 24-position scanivalves. Each scani valve 
contained a 103.42-kN/m 2 (15-psi) differential Statham, variable resistance, unbonded 
strain gage transducer. These transducers are calibrated against a high accuracy 
standard and, if placed in a temperature-controlled environment, will read within an 
accuracy of 0.1 percent of full scale. The transducers were located inside the model and 
subjected to large temperature excursions. During testing in the Boeing Transonic Wind 
Tunnel (BTWT), temperatures recorded at the scanivalves indicated that the accuracy of 
the readout was 0.75 percent of full-scale capability based on the calibration data. For 
tests in the 9- by 7-ft supersonic leg of the NASA Ames Unitary Wind Tunnel, the 
accuracy of pressure measurements was better than ±0.3 percent, based on the 
maximum temperature measured in the test section. 

During the first test in the BTWT (NAS1-12875), two problems were encountered. For 
the first 149 runs, the data filter for one of the scanivalves was inadvertently set at too 
low of a cutoff frequency. This caused a lag that affected five body pressure 
measurements, which produced a maximum error of approximately 0.684 kN/m 2 (0.1 psi) 
at an angle of attack of 16° and M = 0.95. Table A-4 identifies the specific data affected. 
During the first half (approximately) of the test, the scanivalve that recorded lower 
surface pressures between the hingelines for the sections at 2y/b = 0.09, 0.20, 0.35, and 
0.50 was intermittent at an angle of attack of 16°. This problem was eventually traced 
to an electrical problem in the strut. Rather than sacrifice all of these data, the 
incorrect measurements were replaced by extrapolating the data from angles of attack 
of 12° and 14°. 
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In the second test in the BTWT (NAS1-14962), damage to one of the quick disconnects 
for the wing caused the loss of measurements at a series of orifices (x/c from 0.125 
through 0.600) on the lower surface at 0.80 b/2. Replacement values were obtained by a 
linear-span wise interpolation between the adjacent sections. 

WIND TUNNEL FACILITIES AND CAPABILITIES 

The majority of testing (NASA contracts NAS1-12875 and NASI- 14962) of this model 
was conducted in the Boeing Transonic Wind Tunnel (BTWT). There was also limited 
testing (NASA contract NAS1-14141) of the flat and twisted wings in the 9- by 7-ft 
supersonic leg of the NASA Ames Unitary Wind Tunnel. These two facilities are 
described in some detail. 

BOEING TRANSONIC WIND TUNNEL (BTWT) 

The BTWT is a continuous-flow, closed-circuit, single-return facility with an operating 
range of Mach number from 0.0 to nearly 1.1. The test section is 2.438 by 3.658 by 
4.420 m (8 by 12 by 14.5 ft) with 11.0 percent of the wall area in slots. 

The tunnel layout is shown in figure A-7. The tunnel stagnation pressure is 
atmospheric with a total temperature range of 300 K to 356 K (540° to 640° R). The 
variation with Mach number of Reynolds number based on the mean aerodynamic chord 
(M.A.C.) of this model is shown in figure A-8., which also shows the variation of 
dynamic pressure with Mach number. The 26 856-kW (36 000-hp) wound-rotor induction 
motor in tandem with a 13 428-kW (18 000-hp) synchronous motor provides the power to 
drive a 7.315-m (24-ft) diameter fan up to a maximum speed of 470 rpm. The fan is 
made up of a 5.486-m (18-ft) diameter hub with 72 fixed-pitch fiberglass blades 0.914 m 
(36 in.) long in two stages and directs circuit air through two stages of 67 hollow steel 
stators. 

Data System 

The BTWT data system provides the capabilities of real-time test data acquisition, 
feedback control computation, and display. The data system consists of an Astrodata 
acquisition subsystem and a computing subsystem that uses a Xerox data system 
(XDS 9300) digital computer. The Astrodata system acquires signals from the sensors, 
conditions them, and passes them directly to the computer. Test data (averaged from as 
many as 256 samples per test point) are recorded on a rapid-access data drum. As final 
computations are performed, selected on-line displays are provided on analog X-Y 
plotters and tel etype writers. Real-time computations and displays are performed every 
200 milliseconds for control and test monitoring functions. Any test data may be 
retrieved from rapid-access drum storage and displayed on an oscilloscope. On-line 
programs also provide for the preparation of magnetic tapes for plotting or interfacing 
with off-line programs. Figure A-9 is a schematic of the data acquisition and reduction 
system. 

Mach Number 

Mach number in the BTWT is referenced to the horizontal and lateral center of the test 
section at tunnel station 1000, which was the pitch point of this model (40-percent 
M.A.C.). 
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The pressures used in determining the Mach number, p s , and p t are measured through 
permanently positioned sensors. Static pressure p s is measured by a 103.42-kN/m 2 
(15-psi) absolute transducer. A 103.42-kN/m 2 (15-psi) differential transducer is used to 
obtain total pressure by measuring (p t -p s ). These transducers are temperature 
compensated in addition to being in a ±1.11° C (±2°F) environment. Transducer 
performance is checked periodically, and both the static and differential transducers 
have shown a maximum deviation of ±0.02 percent of full scale. 

The static pressure tap is located out of the test section above the ceiling in the 
pressure cap plenum. A correction is made to adjust this static pressure reading to the 
measured test-section-centerline static pressure determined during calibrations at 
station 1000. The tunnel total pressure is obtained from a total pressure probe mounted 
near the tunnel ceiling in the bellmouth throat (fig. A-7). 

Signals from the pressure sensors are fed to the XDS 9300 computer. The XDS system 
computes and updates the Mach display five times per second. Accounting for the entire 
system, calculated Mach number is accurate within ±0.002. Data are recorded only 
when the tunnel is within a preselected Mach tolerance. For this test, a tolerance of 
±0.003 was used. 

Dynamic Pressure 

The dynamic pressure q is computed from the Mach number and the corrected static 
pressure. The estimated tolerance on dynamic pressure is ±95.8 N/m 2 (±2.0 psf). 

Angle of Attack 

The angle of attack of the reference point (0.25 M.A.C. for this model) for a 
sting-mounted model is a combination of the input angle measured at the base of the 
sting and several incremental corrections. The input angle of attack is determined by an 
encoder mounted in the strut. This angle is accurate within ±0.02°. This angle is then 
modified by the effects of sting deflection, up-flow, and wall corrections. 

Sting deflections due to load were determined during the calibrations of the strain 
gages, which are mounted on the integral sting body of the model. These deflections are 
known within ±0.02°. The corrections for sting deflection are based on the normal force 
and pitching moment loads obtained during wind-on data acquisition. The sting 
deflection was taken into account when setting test angles of attack, to minimize the 
variation in final angle of attack for the various model configurations. The strain gages 
attached to the sting body of this model have an estimated accuracy of ±5 percent of 
full-scale reading. This means that the sting deflections based on maximum model loads 
were known within ±0.11°. 

During run 55 of the second test in the BTWT, the wiring for the pitching moment gage 
broke, affecting both the normal force and pitching moment measurements and, 
therefore, the calculation of sting deflection under load. For the remainder of that test, 
the model angle was set using the angle of attack as determined by the encoder for the 
most similar previously run configuration. After the test, the normal force and pitching 
moment obtained by integrating the pressure data were used to correct the final angle 
of attack. To verify this procedure, a comparison of these two methods was made using 
data obtained prior to run 55; the results matched within ±0.01°. 



Up-flow corrections were made based on data obtained from upright and inverted runs 
on a calibration model of similar span. These corrections were less than 0.2°. It is 
generally accepted that the up-flow values are known within ±0.05°. 

A correction to model angle was made for the effect of lift interference for 11-percent 
slotted walls. The lift interference is a function of the ratio of model-to-test section size, 
test section shape, C^, and wall geometry. For Cjj = 1.0, this correction is on the order 
of -0.48°. Due to the limited amount of experimental substantiation, the wall correction 
could be in error by ±20 percent. 

NASA AMES UNITARY WIND TUNNEL, 9- BY 7-FT SUPERSONIC LEG 

The 9- by 7-ft supersonic leg of the NASA Ames Unitary Wind Tunnel is a 
continuous-flow, closed-return, variable-density facility with an operating range of 
Mach number from 1.54 to 2.50. (A schematic is shown in fig. A-10.) The tunnel is 
equipped with an asymmetrical sliding-block nozzle and a flexible upper plate; variation 
of the test section Mach number is achieved by translating, in the streamwise direction, 
the fixed-contour block that forms the floor of the nozzle. For this test, the Reynolds 
number was selected as 8.65 by 10 6 based on the mean aerodynamic chord (c) of this 
model. The test section is 2.74 by 2.13 by 5.49 m (9 by 7 by 18 ft). 

The tunnel air is driven by an 11-stage axial-flow compressor that is powered by four 
variable-speed, wound-rotor induction motors with a combined output of 134 280 kW 
(180 000 hp). Four 850 m 3 (30 000 ft 3 ) spherical storage tanks provide dry air for tunnel 
pressurization. The temperature is controlled by aftercooling. 

Data System 

The data acquisition system is comprised of a Beckman 210 analog-digital recorder and 
a minicomputer. Output from the Beckman 210 is converted to an acceptable format and 
transmitted by the minicomputer to an IBM 360 computer, which is located in the 
AMES Research Center central computer facility for the processing and preparation of 
final data. This flow is illustrated in figure A-ll. 

Angle of Attack 

The angle of attack of the reference point (0.25 M.A.C. for this model) for a 
sting-mounted model is a combination of the input angle at the base of the sting and an 
increment due to sting deflection. The input angle of attack at the base of the sting is 
accurate within 0.02°. 

Sting deflections due to load were determined during the calibration of the strain gages 
mounted on the integral sting body of the model. The corrections for sting deflection are 
based on the normal force and pitching moment loads obtained during wind-on data 
acquisition. The sting deflection was taken into account when setting test angles of 
attack to minimize the variation in final angle of attack for the various model 
configurations. Only a crude calibration of the normal force and pitching moment gages 
was obtained since the force and pitching moment measurements were used primarily 
for calculating sting deflection. Comparison with the integrated pressure results 
indicates that both force and moment measurements may be about 10 percent low, 
which could yield a maximum error in final angle of attack of 0.1°. 
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TESTS AND DATA ACQUISITION 
BOEING TRANSONIC WIND TUNNEL (BTWT) 

Tests 

As previously stated, tests were conducted in the BTWT under two NASA contracts. 
Table A-5 shows the 54 configurations that were tested under contract NAS1-12875. 
The 12 configurations tested under contract NASI- 14962 are shown in table A-6. Two of 
the configurations were included in both tests to ensure that data from the two tests are 
consistent. Photographs of some of the configurations are shown in figures A- 12 
through A-15; a diagram of the model installation in the BTWT is shown in figure A-16. 

Pressure and total force data were obtained at Mach numbers of 0.40, 0.85, 0.95, and 
1.05 for all configurations and at Mach numbers of 0.70, 1.00, and 1.11 for selected 
configurations. Test angles of attack were from -8° to +16° in 2° increments, however, 
not all angles of attack were included for all configurations and/or Mach numbers. 
Tables A-5 and A-6 show the run numbers for each Mach number and configuration for 
which these data were obtained. 

During the first test, wingtip deflection pictures were taken for representative 
configurations at three Mach numbers to evaluate the stiffness of the wing. These were 
compared to wind-off reference pictures to determine the relative deflection and twist. 
Configurations included the flat and twisted wings, and trailing-edge control surfaces 
deflected +30.2°, 0.0°, and - 17.7°. Whereas the tip did deflect (less than 2 cm), the 
change in incidence was negligible even at M = 1.05; and, due to model flexibility, no 
corrections to the data were required. 

Data Repeatability 

Comparisons of data from both tests (NASA CR-165701), show that the data are within 
the tolerances expected for repeat runs during a single test. Therefore, data from both 
tests are compared without regard to the test in which the data were obtained. 

Data Acquisition and Initial Processing 

The pressure data were recorded through the use of fifteen 24-position scanivalves 
located in the fore body of the model. Pressure transducers in the scanivalves measured 
the differential pressure between the local surface pressures and tunnel total pressure. 
Signals from the scanivalves, force and moment data, tunnel parameters, and model 
attitude angle were recorded on the Astrodata system and reduced using the XDS 9300 
computer 

NASA AMES UNITARY WIND TUNNEL, 9- BY 7-FT SUPERSONIC LEG 
Tests 

Table A-7 lists the 13 configurations that were tested. Photographs of two of these are 
shown in figures A-17 and A-18; a diagram of the model installation in the test section 
is shown in figure A-19. Pressure and total force data were obtained at Mach numbers 
of 1.70, 2.10 and 2.50 for all configurations. Table A-7 shows the run numbers for each 
Mach number and configuration for which these data were obtained. Test angles of 
attack were from -8° to +14° in 2° increments and +15°. 
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Table A-5. —( Continued) 










































Table A- 5. —( Concluded) 



423 418 412 407 446 431 438 











































Table A-6.— Summary of Subsonic/Transonic Test Conditions by Run Number 
(NASA Contract NASI -14962) 


■ 

Mach 

number 

Trailing-edge deflection, degrees 

Full span 

Outboard (inbd=0.0) 

Inboard (outbd=0.0) 

0.0 

8.3 

8.3 

8.3 

Twisted wing 

0.40 

15 

25 

30 

20 

0.70 

14 

24 

29 

19 

0.85 

13 

23 

28 

18 

0.95 

12 

22 

27 

17 

1.00 

11 




1.05 

10 

21 

26 

16 

Cambered-twisted wing, fin off 

0.40 

43 

65 

80 

57 

0.70 

41 




0.85 

40 

62 

78 

59 

0.95 

39 

64 

79 

58 

1.00 

38 




1.05 

37 

61 



77 

55 

Cambered-twisted wing, fin on 

0.40 

49 

70 

75 

54 

0.70 

45 




0.85 

48 

68 

73 

52 

0.95 

47 

69 

74 

53 

1.00 

46 




1.05 

44 

67 

72 

51 
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Data Acquisition and Initial Processing 


The pressure data were recorded through the use of fifteen 24-position scanivalves 
located in the fore body of the model. Pressure transducers in the scanivalves measured 
the differential pressure between the local surface pressures and a known reference 
pressure. Signals from the scanivalves, force and moment data, tunnel parameters, and 
model attitude angle were recorded on the Beckman 210 analog-digital recorder and 
reduced by the Ames staff. 

TRIP STRIP 

A trip strip of no. 60 carborundum grit was used throughout the tests with the 
exception of one series. On the body, the trip strip was 0.32 cm (0.125 in.) wide and 
placed 2.54 cm (1 in.) from the nose. On the wing, it was 0.32 cm (0.125 in.) wide from 
the side of body to the midspan control surface break (0.57 b/2), and tapered to 0.16 cm 
(0.0625 in.) wide at the wingtip. On the upper surface of the wing, the trip strip was 
placed at 15-percent chord; and, on the lower surface, it was placed just aft of the 
location of the leading-edge control surface brackets on the flat wing (see fig. A-6). 
Density of the grit was 4 to 5 grains per quarter-inch (6 to 8 grains per cm) of trip strip 
length. 

FINAL DATA 

Final data (pressure coefficients, tunnel parameters, and model attitude) were merged 
on magnetic tapes with appropriate configuration and test point identification for 
integration and plotting of these data. 

A detailed description of the data editing and integration procedure are included in 
appendix B. 




Figure A-2.—Spanwise Twist Distribution for the Model Wing 




Fraction of semi span, 2y/b 

(a) Definition of k, Factor on Basic Camber 


Wing reference 
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Midspan, 2y/b = 0.50 


Wing reference 
plane 




Root, 2y/b = 0.0 


Wing reference 
plane 



(b) Typical Sections 


Figure A -3. — Cambered- T wisted Wing Section Geometry 
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Wing 
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P* anG (a) Undeflected Trailing Edge Control Surface 

Ail dimensions in centimeters (inches) 



Wing reference plane 


(b) Trailing Edge Control Surface Deflected 8.3° 


Figure A-4.-Fin Geometry, Section at 0.725 Semispan 
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plane 
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Figure A-6. -Control Surface Bracket Details 




(b) Test Section 

Figure A-7.— Boeing Transonic Wind Tunnel 
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Figure A-8.— Variation of Reynolds Number and Dynamic Pressure With Mach Number 





Figure A-9.— Data Acquisition and Reduction System— Boeing Transonic Wind Tunnel 
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Figure A- 10.— Schematic of 9- by 7- ft Supersonic Leg of NASA Ames Unitary Wind Tunnel 






Figure A- 7 l.—Data Acquisition and Reduction System— 9- by 7-ft Supersonic Leg 
of NASA Ames Unitary Wind Tunnel 
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Figure A-1 2.— Mode/ in Boeing Transonic Wind Tunnel— Fiat Wing; L.E. Deflection , Full 

Span = 0.0 T.E. Deflection, Full Span = 0.0 ° (NASA Contract NASI -12875) 
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Figure A- 1 3.— Model in Boeing Transonic Wind Tunnel— Twisted Wing; T.E. Deflection Full 
Span = 0.0°; (NASA Contract NASI -14962) 
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Figure A- 14.— Model in Boeing Transonic Wind Tunnel— Cam bered-Twisted Wing, Fin Off; T.E. 
Deflection, Full Span = 0.0°; (NASA Contract NASI -14962) 


56 




Figure A- 15.— Model in Boeing Transonic Wind Tunnel— Cambered-Twisted Wing, Fin On; 
T.E. Deflection, Full Span = 0.0° ( NASA Contract NASI -14962) 
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Tunnel ceiling 
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Figure A-1 6.— Model Installation in Boeing Transonic Wind Tunnel (NASA Contracts NASI -12875 and NASI -14962) 






Figure A-17.—Modei in 9- by 7-ft Supersonic Leg of NASA Ames Unitary Wind Tunnel-Flat 
Wing , Rounded L.E. (NASA Contract NASI -141 41) 



Figure A- 18. —Mode! in 9- by 7-ft Supersonic Leg of NASA Ames Unitary Wind Tunnel— 
Twisted Wing (NASA Contract NASI- 14 14 1) 
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Tunnel ceiling 
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Tunnel (NASA Contract NASI -141 41) 




APPENDIX B 


DATA REDUCTION 


DATA EDITING 

There were some cases encountered with these data where the methods of data editing 
available within the integration programs were not adequate. Because the plotting 
program assumes that the geometry is the same for all configurations, and the 
chord wise location of orifices on the various model parts was not absolutely identical, 
points were added as required. Therefore, some interpolations or extrapolations using 
selected orifices were done before the integration program was used. The row of orifices 
on the body at the wing-body intersection was extended in front of the wing and aft of 
the wing by interpolating between the orifices located at 90° and 135°. 

Some specific problems with the data acquisition systems required the replacement of 
some data for parts of tests. These are identified in appendix A. 

Several methods were introduced into the integration program to replace or add data 
points to account for: 

- Plugged or leaking orifices, or bad data points 

- Extrapolating the data to leading and trailing edges 

- Hingeline discontinuities in the pressure data 

These procedures were selected by code for each point. The codes are described in the 
following list and are illustrated in figure B-l. An additional use of these codes is to 
ensure that only measured pressure data (CODE* = 0) are identified with symbols on 
the plots. The subscript i identifies the position of the point from the leading edge of the 
upper or lower surface of the section. 

IF CODEj = 0, use pressure as entered on tape (measured pressure) 

= 20, use as entered on tape (previously replaced value) 

= 1, interpolate from adjacent points 
= 2, extrapolate from two preceding points 
= 3, extrapolate from two following points 
= 4, set equal to preceding point 
= 5, set equal to following point 
= 6, interpolate using points (i-2) and (i + 1) 

= 7, interpolate using points (i-1) and (i+2) 
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IF CODEj = negative of above, evaluate as above but average with corresponding point 
on opposite surface; used for leading and trailing edges of section only 

Editing of the pressure data is done in the following order: 

1. Each section is done separately. 

2. Each surface (upper or lower) per section is done in the following sequence: 

a. Starting at leading edge, points with codes of 1, 2, and 4. 

b. Starting at trailing edge, points with codes of 3, 5, 6, and 7. 

3. Leading- and trailing-edge points with negative codes are evaluated. Both upper 
and lower surface codes need not be negative and need not be the same negative 
code. 


CALCULATION OF NET PRESSURE COEFFICIENTS 
The net lift distribution on the section is calculated by: 


r = f 

p,net ^p, lower ^p, upper 


(B-l) 


INTEGRATION OF PRESSURE DATA 


To account for the effects on integrated coefficients of the deflected control surfaces, 
each streamwise section (of which there are NSECT) is divided into segments (of which 
there are NSEG). These segments are the leading-edge control surface, wing box, and 
trailing-edge control surface. The upper and lower surfaces of each are integrated 
separately over the number of points available ((number of orifices + 2) = NP1), and 
are based on the segment chord length c. Sign conventions are shown in the following 
sketch. The equations, which use a rectangular integration process, follow. 
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SEGMENT COEFFICIENTS 

Integration of the pressures for each segment per surface per section is the first step. 


o Normal force coefficient C 


n,s 


NP1 


Cn,s = °5 E [( C p), + ( C p) M ][©r(r|,] 

C n,s,net = ^n,s,lower “ C n,s, upper 


o Chord force coefficient C 


c,s 


C c,s *0.5 E[(C p ) j+( C p)| _ ,][(§>-(!>_] 


^c,s,net ^c,s, upper " ^c,s,lower 

Pitching moment coefficient about segment leading edge C m £ 

NP1 


C m,s = °- 5 E [( C p)i + ( C p) w ' 


-* + 
l_ C i-1 




2.0 




NP1 


= (U5 g [( C p), + ( c p)m’ 


(f) -(f) 


i-1 


C m ,s,net " ^m,s, upper " C ms> i ower 


o Pitching moment coefficient about 0.25 c of segment C m 25c, s 


^m.25c,s ^m,s, + ^-25 ^n,s 


(B-2) 

(B-3) 

(B-4) 

(B-5) 

(B-6) 

(B-7) 

(B-8) 
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SECTION COEFFICIENTS 

Total section coefficients are obtained by summing the segment coefficients, taking into 
account both the segment deflections as defined in the following sketch and the segment 
chord lengths. These coefficients are based on the section chord length c. 


-8 



This sign convention for leading-edge deflection is used only in summation 
of coefficients. 



o Normal force coefficient C n 


NSEG . , NSEC , . 

o Pitching moment coefficient about section leading edge C m 


(B-9) 


m 


NSEG / \2 r l/c\ 

= Z (< m . S ). (i) + ( C n,s) , (I- 0 " 6 1) + ( C c,s) , si " «1 (i) 

j— 1 J ' 'i L J 


NSEG / \ rv _ v 

- S [( c n.s) j «* *j - (Cc.*)j 6 j] (i\ [ L E c' L E j ] j 

J - J 


(B-10) 


where 

c s 

c 

8 


x L.E.,s 

X L.E. 


is segment chord length, cm 
is section chord length, cm 

is deflection of segment relative to section chord plane, leading edge 
up, degrees 

is leading edge of segment, cm 
is leading edge of section, cm 


o Pitching moment coefficient about 0.25 c of section C m 2 s c 

C m .25c = C m +0 - 25 C n 


(B-ll) 
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TOTAL SURFACE COEFFICIENTS 

To obtain total surface coefficients, the assumption is made that the section coefficients 
apply for a finite distance on both sides of each row of orifices. The equations for total 
surface coefficients are as follows: 

o Normal force coefficient Cjj 

NSECT 

Cn= ^ ( Cn )k( Sh )k (B-12) 

o Bending moment coefficient Cg 

NSECT 

C B~S(b/2) ^ ( Cn )k ( S b y )k (B-13) 

o Pitching moment coefficient about 0.25 M.A.C. C^j 


NSECT / 

c M=^ S |( C m) k ( S h c ) k + ( C n) k ( S h) k [ x ref-( x L.E.) k ] 
where 


c is reference chord for pitching moment, cm 

x ref is reference station for pitching moment, cm (0.25 M.A.C.) 

Xl.e. is leading edge of section chord, cm 

b/2 is reference length for bending moment, cm 


(B-14) 
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DETERMINATION OF GEOMETRIC CONSTANTS REQUIRED FOR 
INTEGRATION 

To obtain total surface coefficients, the assumption is made that the section coefficients 
apply for a finite distance on both sides of each row of orifices. The input geometry 
required to calculate the areas, and products of area and length required for the 
summation of total surface coefficients, is shown in the following sketch. 



o Section area: 


f Vout r x t + tan A T.E.y 
J y-m •' x L + ,an A L.E.y 


dy dx 


= c r( y out " y in) + 0 5 ( tatl A T.E. " tan A L.E.) (yout 2 " Vin 2 ) (B-15) 


o Product of section area and mean chord: 

.y ou t r x T + tan A T.E.y 


s h c = 


/-yout r 

I c dy dx 

y; n ' /x L + tan A L.E.y 


= c r 2 (y out - yj n ) + c r ( tan A T.E. ' tan A L.E.) (y 0 ut 2 ' ^ 11 ”) 


(B- 16 ) 
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o 


Product of section area and bending moment arm: 


/•Vout r X T + tan A T.E. y 
s h y = / / (y - Vreljdy dx 

"yin Jy ‘L + tan A L.E.y 

_ c r -(tan A T E .-tan A L E ,)y re f / 2 „ 2\ 

- 2^0 \ y out y in / 

+ tan A T E ^'J an A L . E . (yout 3 . ^3) . Cf yref ( Yout - y J 
o Product of section area and leading-edge coordinate: 


/•y out /• A T T iajl yv " 
•'vin *x L + tan A UE . 


x T + tan Ay £ y 


ShX 


x l.e. d y dx 


/ \ tan A L.E. c r + x L( tanA T.E.- ,an A L.E.V. 2 .. 

= x L c rV y out-yin^ + lo V°ut ' Vi 


in 


+ tan A 


(tan A T E - tan A l .e) /-.. 3 v 3\ 

L.E. 3.0 \ y out y in / 


o Total surface reference area: 


NSECT 

s = L (S h ) 

k= 1 v n 


M.A.C. and X coordinate of M.A.C. leading edge: 

NSECT 

_ 1 


c = s E (v) 

5 k=l 7 


NSECT 

X L.E.. M.A.C. =4- E ( S h x )i 


k= 1 


The required integration constants for the wing and body are shown in table B-l. 


(B-17) 

0 

(B-18) 

(B-l 9) 

(B-20) 

(B-21) 
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Table B-1.— Integration Constants 


Reference area = 3128.45 cm 2 
M.A.C. = 75.311cm 

Half span = 50.80 cm 


Pitching moment referenced to 0.25 M.A.C. 

Bending moment referenced to 0.086 ~ (y ref = 4.374 cm) 


Wing 


L.E. of M.A.C. at B.S. 87.760 cm 


2y/b 

Ay 

(b/2) 

Area 

cm 2 

Area . chord 
cm 3 

Area • (y-V re f) 
cm 3 

0.09 

0.0425 

219.69 

22 357. 

167. 

0.20 

0.1575 

733.51 

67 415. 

4 206. 

0.35 

0.1500 

580.54 

44 374. 

7 857. 

0.50 

0.1400 

437.93 

27 084. 

9 t48. 

0.65 

0.1600 

377.64 

17 722. 

10 729. 

0.80 

0.1300 

210.35 

6 794. 

7 528. 

0.93 

0.1400 

129.79 

2 487. 



5 505. 


Body 


Longitudinal 

section 

Area 

cm 2 

Area • L 

cm 3 

1 

356.61 

81 258. 

2 

504.32 

114916. 

3 

70.94 

16 164. 
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Extrapolated 
using code = 2 




When 

code = —2 
(similar 
at leading 
edge when 
code = —3) 


Figure B-1.— Codes Used to Interpolate and Extrapolate 
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APPENDIX C 


METHOD FOR CALCULATING NET LEADING-EDGE SUCTION 

COEFFICIENTS 

The procedure is described for obtaining leading-edge suction coefficients (and therefore 
leading-edge vortex strength) from the net pressure distribution that i3 predicted by the 
linear theory FLEXSTAB. Also explained is a procedure to account for the retarding 
effect of a rounded-leading edge. 

As the pressures for a specific wing shape are used, the effect of wing warp and twist 
are automatically included. 

GENERAL LEADING-EDGE SUCTION CALCULATION PROCEDURE 

The leading-edge suction coefficient C B at each spanwise location can be calculated from 
the local wing geometry and net lifting pressure distribution as: 


C 


s 


2.0 7T 

be 


cos A 

[2.0 S ref J 



tan A k C 


2 

l 


where 

A 

b 

c 

S ref 

k 

M 

Ci 


= the local leading-edge sweep angle 
= wing span 
= local chord length 
= reference wing area 



= Mach number 


the strength of the net lifting pressure coefficient 
singularity near the leading edge of the wing 

SS %r ir 


Cp.net = ne t lifting pressure; wing lower surface pressure 
minus upper surface pressure 

i = (x - x LE )/c, fraction of local chord 

x - x L e = local streamwise distance from the wing leading edge 


(C-l) 


(C-2) 
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The net lifting pressure coefficient can be approximated by a simple fractional 
polynominal equation in terms of the fraction of local chord £ as: 


^p^et 

4.0 


£l_ 

VT 


+ c 2 T /2 + C 3 i^/2 + 


(C-3) 


This equation provides a direct relationship between the net lifting pressures and the 
leading-edge pressure singularity, C x . 

In principle, Cj could be determined by a least squares curve fit of equation (C-3) to the 
C p> net distribution calculated by the FLEXSTAB program. However, the 
FLEXSTAB-calculated pressures near the wing leading edge are uncertain because of 
the inability of a constant-pressure-panel method to properly represent the leading-edge 
singularity. The pressures do tend to converge to the exact values farther aft on the 
wing chord. The irregularities of the calculated pressures make it difficult to accurately 
determine Cj by this procedure. The fact that integration tends to smooth numerical 
irregularities of a calculated function, is utilized to calculate C x . 

Integrating equation (C-3) provides an approximate function for the local integrated 
load distribution as: 


/ 


-'p.net 

"4ir 


df = 2.0 Cj vT + f§ c 2 £ 


3/2 


+ || C 3 


+ 


(C-4) 


Define the function g(C p ) as: 


g(Cp) 



°p,net 

4.0 


d^ 


2.0 vT 



o 


8.0 VT 


(C-5) 


Combining equations (C-4) and (C-5) provides an approximating function for g(C p ) as: 


g(C p ) = Ci+|| c 2 f + i|c 3 


5.0 




(C-6) 


The leading-edge pressure singularity Cj can be determined by integrating the 
FLEXSTAB pressures successively back to th^ aft edge of each panel to obtain the 


71 


distribution of g(C p ) over the wing chord. The approximating equation (C-6) is then "fit” 
to the calculated values of g(C p ) by a least squares approach. 

Typically, if the values of g(C p ) back to approximately 25-percent chord are used, a 
simple quadratic equation for g(C p ) is adequate (i.e., use the first three terms on the 
right-hand side of equation (C-6». 

FLEXSTAB uses constant pressure panels to represent the wing pressure distribution. 
Therefore, the value of g(C p ) back to the aft edge of panel j (that is, for |j) can be 
represented by a summation: 


/ f J 1 

^p,net [pp.net] j j^i ~ £i-lj 


g(C D ) = 


8.0 


8.0 


(C-7) 


The procedure for this evaluation is illustrated in figure C-l. It should be noted that the 
fraction of local chord £ used in this evaluation is that of the aft edge of each panel. 

The leading-edge suction distribution for a specific combination of wing shape, Mach 
number, and angle of attack can be calculated as follows: 

1. Calculate the net lifting pressures for the specified wing geometry at the desired 
Mach number and angle of attack. 

2. Calculate the corresponding g(C p ). values at each semispan station where the 
pressures are obtained, using equation (C-5) or (C-7). 

3. "Least squares fit” equation (C-6) to the calculated g(C p ) values on the forward 25 
to 30 percent of each wing section to determine C-l tat least five or six values of 
g(C p ) should be used). 

4. Use equation (C-l) to determine the local leading-edge suction coefficient. 

ANGLE OF ATTACK SCALING 

FLEXSTAB lifting pressures vary linearly with angle of attack. The local leading-edge 
singularity at any angle of attack can be obtained by linear interpolation from any two 
known data sets for the same wing shape and Mach number. The full leading-edge 
suction coefficient can then be calculated using equation (C-l). 
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LEADING-EDGE VORTEX DEVELOPMENT AND NET SUCTION 
CALCULATION 


A method for predicting the retarding effect of rounded-nose airfoils on the progressive 
development of leading-edge vortices on highly swept wings is presented in references 
18 and 19. The leading-edge vortex is assumed to be present at any station on the wing 
when the local leading-edge suction coefficient C 8 is greater than the local airfoil nose 
pressure force coefficient. 

It is shown in reference 19, that the leading-edge vortex starts at a particular spanwise 
station t? when 

Ci = C 1>R 


where: 


C 


1,R 


Rl 


~ V2.0 c 
k tan A 


(C-8) 


is the ratio of the airfoil nose radius (perpendicular to the wing leading 
edge), to the local streamwise chord. 


The perpendicular nose radius to streamwise chord ratio can be obtained approximately 
from the local streamwise chord ratio R s /c as 


R-L _ 1,0 Rs 

c cos A c 


(C-9) 


The net leading-edge suction coefficient is calculated using equation (C-l) by replacing 
the term C x with the term Cj - Cj R . 
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Panel 

* 

Cp,net 

g(c p ) 

1 

.025 

1.7942 

.035461 

2 

.050 

.7461 

.035502 

3 

.075 

.6064 

.035907 

4 

.100 

.5126 

.036162 

5 

.150 

.4507 

.036799 

6 

.200 

.3865 

.037270 

7 

.250 

.3497 

.037707 

8 

.350 

.3204 

.038638 

9 

.450 

.2858 

.039401 

10 

.550 

.2619 

.040054 

11 

.650 

.2398 

.040562 

12 

.700 

.2180 

.040715 

13 

.750 

.2066 

.040826 

14 

.800 

.1897 

.040855 

15 

.900 

.1656 

.040700 

16 

1.0 

.1047 

.039920 


.04 

- 

T 1 1 — 1 

0 o__- 

r"°' \ 

.03 

C 

\ Calculated 

\ (eq.C-7) 

g(c p ) 

.02 

- 

Curve fit, first 7 points 

.01 

- 

- 

0 

1 — 

j 1 1 1 1 1 1 1 


0 .2 .4 .6 .8 1.0 


Fraction of local chord, £ 


A least-squares curve fit of 
equation (C- 6 ) yields: 

Ci = 0.035059 
C 2 = 0.035244 
C 3 = - 0.020560 


Figure C- 1.— Sample Calculation of C /, Flat Wing ; M = 0.85, a = 10°, 7 ? = 0.3 19 
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Figure 2.-Spanwise Twist Distribution for the Model Wing 


Wing reference 
plane 


Tip, 2y/b = 1.00 



Midspan, 2y/b = 0.50 


Wing reference 
plane 


Root 2y/b = 0.0 


Wing reference 
plane 



Figure 3.~Cambered-Twisted Wing Section Geometry 
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Figure 4. -Pressure Orifice Locations 





Figure 5. -Upper Surface Isobars , Flat Wing, Rounded Leading Edge, M = 0.85 


Section normal force 
coefficient, C n 



Flat wing 

Cambered-twisted wing 


Figure 6. -Comparison of Experimental Section Normal Force Coefficients for Two 
Wing Shapes, M = 0.85, 77 =0.80 





Error 


O Flat wing experiment 
A Cambered-twisted wing experiment 
Predicted cambered-twisted wing 


(a) Prediction Using Linear Theory, FLEXSTAB 



(b) Prediction Using Linear Theory, FLEXSTAB, for Increment to Flat Wing Data 


Figure 7.— Typical Procedures Used for Aerolastic Loads Prediction, M = 0.85, r? = 0.80 
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Angle of attack, a, degrees 


Fraction of semispan, 17 = y/(b/2) 


(a) Flat Wing 



0.4 



Fraction of semispan, 17 = y/(b/2) 


(b) Twisted Wing 


Figure 8.— 


Net Leading-Edge Suction Coefficients for the Three Wing Shapes, Rounded 
Leading Edge, M = 0.85, Calculated From FL EXSTAB Pressures 












Establish J 
data base \ 


Iterative 

nonlinear 

aeroelastic 

solution 


v 



Figure 9. —Proposed Iterative Aeroelastic Loads Prediction Method Based on Section 
Normal Force Coefficient 
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Net pressure coefficient, C p Met pressure coefficient, 


77 = 0.09 

O Twisted wing 
A Flat wing 

C n = .435 

^ 





















Net pressure coefficient, Cn Net pressure coefficient. 



0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 

Fraction of local chord — x/c 


Figure 1 l.—Chordwise Pressure Distributions Which Produce the Same Section Normal 
Force Coefficient, Twisted and Cambered-Twisted Wings, M = 0.85 
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Net pressure coefficient, Cn Net pressure coefficient. 




Figure 1 1.— (Continued) 
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.8 


.6 - 

.4 - 

Section normal force 
coefficient — C n 

0 - 
-.2 - 
-.4 - 

-8 -4 0 4 8 12 16 

Angle of attack, degrees 

.6 
A 
.2 

Section normal force 0 
coefficient — Cn 

-.2 
-.4 
-.6 
-.8 

-8-4 0 4 8 12 16 

Angle of attack, degrees 

.4 
.2 
0 

Section normal force ^ 
coefficient — C n 

-.4 

-.6 

Wing 

O Flat g 

□ Twisted 
Camber-twist 

Figure 12.— Comparison of Experimental Section Normal Force Coefficients for 
Three Wing Shapes, M = 0.85 
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Figure 14.— Definition of Crossflow Wing Sections 
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Rigid airplane 
model geometry 


Wind tunnel 
pressure data 


Linear 

aeroelastic 

solution 


c p,NL vs a 


I 1 

>4 Linearize I- 1 

1 1 


Cp.lin. vs ot 


^Cp,NL vs a 


Data base 


^p,NL vs C Sf net 


Nonlinear 

aerodynamic 

analysis 


Airplane 

geometry 

Flight 

condition 


FLEXSTAB 
aeroelastic analysis 
Cp = Cpjm. + ^p,NL 


Cs,net 


Structural design loads 


Figure 15.— Proposed Iterative Aeroelastic Loads Prediction Method Based on Net 
Leading-Edge Suction Coefficient 




^Cp,NL " Cp,exp. ~ Cp,lin. 


Figure 16.— Procedure to Establish a Data Base Which Relates Nonlinear Pressure 
Coefficients to Net Leading-Edge Suction Coefficient 
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a; degrees a; degrees 

(a) r? = 0.20 



o; degrees a; degrees 

(b) 7?=0.35 


Figure 18.— Prediction of Streamwise Section Coefficients for the Cambered-Twisted 
Wing Using the Data Base for the Flat Wing Data Related to the Net 
Leading-Edge Suction Coefficient , Compared to Linear Theory, M - 0.85 
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(f ) 77 = 0.93 


Figure 18.— (Concluded) 
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Section pitching moment coefficient — C m .25c Section normal force coefficient 


.8 



T] = .93 


\77 = .35 


Figure 20.— Prediction of Stream wise Section Coefficients for the Cambered-Twisted Wing 

Using the Data Base for the Twisted Wing Data Related to the Net Leading-Edge 
Suction Coefficient, M = 0.85 


I 
















O Flat wing 
A Twisted wing 
□ Cambered-twisted wing 


x/c = 0.450 



Figure 21.— Incremental Nonlinear Net Pressure Coefficients as a Function of the 
Net Leading-Edge Suction Coefficient for the Three Wings, r\ = 0.50, 
M = 0.85 
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Section moment coefficient — C m Section normal force coefficient 



O Flat wing 
□ Twisted wing 
^Cambered-twisted wing 



Note: 

Constant body station (crossflow} section; 
intersects leading-edge at 77 = 0.957 


Figure 22. — Typical Variation of Crossflow Section Coefficients with Angle of 
Attack for the Three Wings, M = 0.85 
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77 = .519 


.8 



(a) Section Normal Force Coefficient 

Figure 23. — Prediction of Crossflow Section Coefficients for the Twisted Wing Using 
the Data Base for the Flat Wing Data Related to the Net Leading-Edge 
Suction Coefficient ; M = 0.85 
























.9 


.9 



\/Cs,net 



\/Cs,net 





net 


O Experiment 
A Predicted 


(b) Section Moment Coefficient at Side of Body 

Figure 24.— (Concluded) 


109 













(a) Section Normal Force Coefficient 


O Experiment 
A Predicted 


Figure 25.— Prediction of Crossflow Section Coefficients for the Cambered-Twisted 
Wing Using the Data Base for the Flat Wing Data Related to the Net 
Leading-Edge Suction Coefficient, M = 0.85 
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(b) 77= 0.319 


Figure 26.— Prediction of Crossflow Section Coefficients for the Cambered-Twisted 
Wing Using the Data Base for the Flat Wing Data Related to the Net 
Leading-Edge Suction Coefficient, Compared to Linear Theory, M = 0.85 





Cn 


Cn 



a, degrees 




Cm 



G; degrees 


(d) 77 ® 0.519 


Figure 26.~(Continued) 
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<e) t?= 0.614 



4 8 12 16 

a, degrees 



(f) 77 = 0.703 


Figure 26. —(Continued) 









(g) 77 = 0.793 




(h) 77= 0.878 




(i) 77 = 0.957 


Figure 26.— (Concluded) 












Net pressure coefficient. 



Figure 27.— Pressure Distribution on the Cambered-Twisted Wing Based on the 

Leading-Edge Suction Coefficients at Constant Body Stations and the 
Flat Wing Data Base , Stream wise Section , 7 7 = 0.50 , M = 0.85 f a= 12° 
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